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^ ■ Abstract 

Q^ , We calculate the complete matrix of two-body scattering amplitudes for the 

s.^ . scattering of longitudinally polarized gauge bosons W^, and Higgs bosons 

^f) . to two loops in the high-energy, heavy-Higgs limit ^/s S> Mh 3> Mw- Use of 

<^ . 

, the Goldstone boson equivalence theorem reduces the problem to one involv- 



' ing only the scalar fields w^, 

and the Higgs boson. Renormalization of the scattering amplitudes requires 



, ing only the scalar fields w^, z (the Goldstone bosons of the original theory) 



Oh! the calculation of the self-energy functions n^(-Mj?), the renormalization con- 

stants and the bare quartic Higgs coupling Aq to two loops. The results 



^ \ will be useful in other calculations. To facilitate the calculations, we introduce 

H : 

, 5^ , a powerful new technique for evaluating integrals over Feynman parameters 

in dimensional regularization which is based on a Barnes' type representation 
of the binomial expansion. We also collect some useful integrals which extend 
the tables given by Devoto and Duke. 
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Typeset using REVT^ 



I. INTRODUCTION 



It was pointed out some time ago by Dicus and Mathur and Lee, Quigg, and Thacker 
that the presently unknown mass Mh of the Higgs boson could be bounded above in 
a weakly-interacting standard model by using the constraint on the magnitude of partial- 
wave scattering amplitudes implied by unitarity. This was examined in detail in where 
it was shown that the only real constraints from two-body scattering arise from processes 
which involve only longitudinally polarized gauge bosons , Zl and Higgs bosons H. 
These processes are enhanced by factors of Mjj/M^^ oc X/g"^ relative to those which involve 
transversely polarized gauge bosons and the small gauge couplings g. Here A is the quartic 
coupling for the Higgs boson, and is related to Mh and the vacuum expectation value 
V = (v^G^)"^^^ = 246 GeV by A = Mjj/2v^. When A (or Mh) is large enough, the 
interactions of the Higgs sector of the standard model become strong, and the unitarity 
constraints on A or Mh are to be interpreted in terms of a transition from weak to strong 
coupling rather than as upper bounds on either. 

The original tree-level analysis of scattering in the neutral channels W^W^ , Z^Zl, HH, 
and ZlH given in has been extended to one loop by a number of authors ||5|-|9|JTT|-|13| , and 
sharpened at tree level P,p!0| and at one loop |jll],|13[ using renormalization-group arguments. 
The most restrictive analysis was that in pT]JT3|] where it was shown that strong coupling 
sets in for A ^ 2.2-2.5, a value substantially below the tree-level bound A ~ Stt/S ~ 8.38 
(or Mh = 1.007 TeV) derived by Dicus and Mathur and Lee, Quigg, and Thacker 0. 
It is desirable to extend the analysis to two loops. This is important as a check that the 



conclusions of |[11|J13|1 are in fact correct, and that an increase of one order in perturbation 
theory does not loosen the (unexpectedly low) upper bound on A in a weakly interacting 
theory. 

The present paper is devoted to the calculation of the two-body scattering amplitudes for 
W^, Zl, H scattering to two loops in the high-energy, heavy-Higgs limit Mw <^ Mh <^ v^- 



Our results are used in the following paper to reexamine the limits on A and Mh to 
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two loops in a weakly interacting theory. 

The organization of this paper is as follows. In Sec. 0, we discuss the general formulation 
of the problem, its reduction through the use of the Goldstone boson equivalence theorem 
to a problem involving only the scalar sector of the unbroken standard model, and the 
renormalization of the scalar Lagrangian. In Sec. fT|, we discuss the calculation of the self- 
energies and wave function renormalization constants to two loops. The most difficult parts 
of the calculation are in the extraction of the finite parts of the two-loop terms. We sketch 
the methods which we found useful for various diagrams, but give only minimal details. 
The results are summarized by diagram in Appendix 0, while some technical points of 
the methods used are discussed in Appendices B, C and D. In particular, we illustrate in 
Appendix B a powerful new method for evaluating integrals over Feynman parameters in 
dimensional regularization which is based on a Barnes'-type representation of the binomial 
expansion. We sketch the use of Kotikov's method of differential equations to evaluate 
certain two-loop graphs in Appendix C, and collect a number of useful integrals in Appendix 
D. 

The complete self-energy counterterms and the wave function and coupling renormal- 
izations are collected in Sec. [llip. These and the renormalized two-loop scattering ampli- 
tudes given in Sec. IV are the principal results of the paper, and constitute the input for 



1^ . The counterterms and renormalization constants are also directly relevant to calcu- 
lations of W^, Zl scattering at "low" energies for very massive Higgs bosons, that is, for 

II. GENERAL FORMULATION OF THE PROBLEM 

A. Framework 

We will be concerned in the following sections with the calculation to two loops (0(A'^)) of 
the complete set of two-body scattering amplitudes for high-energy scattering in the coupled 
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neutral channels W^W£, ZlZl, HH, HZl. The calculation can be simplified enormously in 
the limit we will consider, 3> Mh ^ Mw, by the use of the Goldstone boson equivalence 
theorem [^^|T5|-0. This theorem states that the scattering amplitudes for n longitudinally- 
polarized vector bosons W^,Zl and any number of other external particles are related to the 
corresponding scattering amplitudes for the scalar bosons w^, z to which W^, Zl reduce 
for vanishing gauge couplings g,g' by 



T{Wt Zl, H,...) = {iCTT{w^, z,H,...) + 0{Mw/V^). 



(1) 



The constant C depends on the normalization scheme used in the calculation |]T6|,|T8|,|l9[ , 

1/2 



(J _ ( 



1 + 0{g' 



(2) 



Mw \Z^ 

where the Z's are the wave function renormalization constraints for the physical fields 
and the scalar fields w^. C is equal to unity for in schemes in which the renor- 

malizations are defined at mass scales m <^ Mn. We will renormalize the w'^ , z fields at 
m = 0, a choice which corresponds to massless Goldstone bosons and is consistent with the 
assumption that M^ji/ -C M^. Then [^,0, 



C = 1 + 0(^7^ 



(3) 



Since we are interested in upper limits on A or M^-, the limit M'^/Mfj oc g"^ /\ <^ 1 required 
for the validity of the theorem is natural for our purposes. We will henceforth set the gauge 
couplings to zero, and work entirely with the scalar fields. 

B. Lagrangian and renormalization 



The original Lagrangian for the symmetry-breaking sector of the standard model. 



(4) 



with D^' = 8^" + igW^ ■ T + WB^ and 
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Wi+iW2 ^ 

$ = , (5) 

^ h +iw3 ^ 

reduces for g,g' to an S0(4)-symmetric Lagrangian involving only the scalar fields h 
and w = {wi,W2,ws). The field h will be taken as usual as the component of $ which 
acquires a vacuum expectation value, u at tree level and v in the full theory. We will write 
h as h = V + H where {yL\H\VL) = with respect to the physical vacuum Adding 
all possible S0(4)-symmetric counterterms and rewriting the potential in Cq in terms of v 
rather than u, we obtain the complete Lagrangian for g = g' = 0, 

$ is now to be written in terms of w, H, and v. 

The S0(4) symmetry of the unbroken scalar theory and the condition {Q\H\Q) = 
impose constraints on u"^, A, Sm^, and 6X. For example, Higgs tadpole diagrams shift 
the vacuum expectation value and contribute apparent mass terms for the Goldstone bosons 
which must be absent in the final theory. Taylor has outlined a procedure which sums the 
tadpole diagrams to all orders, and enforces the vanishing of {Q\H\Q). The same condition 
guarantees the vanishing of the masses of the Goldstone bosons w, and allows the tree- level 
vacuum expectation value u to be eliminated in terms of v, the self-energy function 11^(0), 
and the remaining parameters The parameter A appears in the kinetic terms for both 



the w and H fields because of the underlying symmetry in the problem. It is related to the 
wave function renormalization constant hj = 1 + A. Finally, (5m^ and SX are related 
through the constraint Mfj = 2Xv'^, which is taken to be exact, that is, to hold to all orders in 
perturbation theory, with Mh the physical mass of the Higgs boson defined by the condition 
Re (iA]^^(M|f)) = In enforcing this constraint, we are choosing a renormalization 

scheme and definitions of Mh and A which are consistent with those introduced by Sirlin 
and by Sirlin and Zucchini [^] for the full theory, with g,g' ^ 0. The result of this 



analysis ||2T]] is a modified Lagrangian 



£ = Wo ■ S'^wo + Id^H, d'^Ho -\{mI- Re n?,(M^)) Hi + (0)w2 

(7) 



We have introduced bare fields wq = Z\l'^\^ and Hq = Zl/"^!! to keep the kinetic terms and 
free propagators in standard form [Q. 11^ and 11^ are the self-energy functions for the 
bare fields and are to be calculated using the Lagrangian in Eq. (7). We will later drop the 
superscript on 11° for simplicity in denoting the one- and two-loop contributions to these 
functions, 11'' = 11^^) + 11*^^) + ■ ■ ■. The scattering amplitudes calculated in terms of wq and 

1/2 

Hq require the usual wave function renormalization, so must be multiplied by a factor 
for each external line. 

The form of the new Lagrangian clearly ensures that the Goldstone bosons w are massless 
and that the Higgs boson has the physical mass Mh- The constraint Mfj = 2Af^ can be 
used to eliminate either Mfj or A given the physical value of v. However, it will be most 
convenient to keep C as written, and only use the constraint as needed. Finally, in using 
this form of the Lagrangian, all Higgs tadpole graphs should be omitted since the shift in 



the vacuum expectation value from m to w has been suppressed following Taylor . 

The SO (4) symmetry of the original Lagrangian is broken by the mass terms and the 
trilinear couplings in Eq. (|^, but remains in the quartic couplings. There is still a residual 
SO (3) symmetry for the w fields, so the self-energy functions and renormalization constants 
for the charged and neutral particles w"^ = {wi =F ^^2)/^/^ and z = W3 are identical, 

z:' = Z-' = 1 - ^Xip%.=,. (8) 

However, the self-energy function and wave function renormalization constant for the Higgs 
boson are different, and the subtraction point is at = rather than = 0, 

Z^' = 1 - ^ReIlW)\p^=Ml,. (9) 

The SO (4) symmetry of the theory is partially restored in scattering processes in the high 
energy limit y/s ^ Mfj- In particular, scattering graphs which involve the dimension-three 
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trilinear couplings in Eq. (|^) are suppressed by powers of s, t, or u relative to graphs which 
involve only the quartic dimension- four couplings |^5[, and can be neglected for s ^ 
subject to some subtleties discussed later. The only contributions of the trilinear couplings 
for ^/s oo occur through the renormalization constants Z^, Zh which are fixed at low 
energies or, at two loops, in the unitarity sum for the partial-wave amplitudes, through 
2 — ^> 3 processes which have a low-energy component. The Z's and low-energy processes 
break the SO (4) symmetry even though it is present in the high-energy interactions. 



C. The Feynman diagrams 

We used the programs DIAGRAMMAR and DRAW |^ to generate the self-energy and scatter- 
ing diagrams which are needed through two loops, and to calculate their weights, including 
the couplings. The results were checked by hand and by a Wick expansion. The thirty-eight 
distinct diagrams for the one- and two-loop self-energies are shown in Figs. without 
the couplings or weights. All tadpole graphs have been omitted as discussed above. The 
coupling counterterms are not shown as these do not alter the Feynman integrals to be 
calculated. 

For ease of reference, the distinct topologies for the self-energy diagrams are named Scoop 
(Q), Bubble {B), Triangle (T), Acorn {A), Eye-in-the-Sky (S), Sunset {S), and Lemon {C). 
Individual graphs are labeled by the number of internal Higgs lines as shown in the figures, 
with extra labels m, 0, and Q to distinguish otherwise-ambiguous cases. In the analytic 
expressions developed in the following sections, the symbols Q, B, ... will refer specifically 
to the Feynman integrals associated with the corresponding diagrams, with the coupling 
constants and weights from the Wick reduction omitted, e.g., from Fig. |^ 

/i i 
(P-M^)((A; + gP-M^) ^^^^ 

for incoming momentum q. 

The scattering diagrams which are needed through two loops are shown in Fig. ^, with 
a similar labeling. The tree, one- and two-loop graphs V, Bs, T>s and As involve only the 
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S0(4)-symmetric quartic couplings. The complete set of scattering diagrams also includes, 
where appropriate, contributions from the t and u channels. These scattering amplitudes 
are S0(4)- and crossing-symmetric and increase in magnitude for s —>■ oo. In contrast, the 
sample diagrams shown in Fig. ^, which include trilinear couplings, vanish for s, \t\, \u\ oo 
as is easily checked. 

The calculation of the self-energy terms will be discussed in the next section. The 
scattering graphs will be considered in Sec. IV. 

III. CALCULATION OF SELF ENERGIES AND RENORMALIZATION 

CONSTANTS 

A. The general problem 

The self-energy function 11^ (p^) for the Goldstone bosons w is given to two loops by 

_ .^0 _ z!^ A _ Reng(M^)-nW(o) 

i-ikvl ( _ Reng(M^)-nW(0) 
\ Ml 

-i\BQ\iReYi^^\Ml) - 23AQ] - i\{-i\vf[l%AiQ + 6^2q + 20Aq] 
+ {-i\v)H[{iReYi%\Ml) - i3\Q)T2 + (m^(0) - i\Q)Ti] 

-iX{-iXvy[40Aio + 24^3 + IQA20 + 8(i3i)2] 

+{~tXvy[72S^ + 2AS; + 16^2 + 48£3 + le^o], (H) 

where Q, B, etc. refer to the Feynman integrals which are associated with the diagrams in 
Figs. We have not given the explicit expressions for the factors in the first two terms 
which arise from the wave function and coupling renormalizations so as to indicate their 
origin more clearly. These factors can be expanded systematically using the one-loop results 
from Eq. (11) and the following equation and their analogs for dU^/dp"^ and dU^l/dp^. 
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The corresponding result for the Higgs boson is 

'Z'w \ j 



-iX{-iXvf[5AAiQ + I8A2Q + 12Aiq] 

+ {-iXvf[36{iReU'H\Ml) - i3\Q)T^ + 12(zn^(0) - i\Q)%] 
+{-i\f[QSi + 6^3] - i\{-i\vf[2lQAi + 72^2m + 48Am] 
-i\{-i\vf[3Q{BQf + 36i3oS2 + 54(^2)'] 

+ (-zAt;)^[648^5 + 216^3 + 48£i + 648£5 + 144£2m + 24£i]. (12) 

We will only need these functions evaluated for = (11° (0)) and = Mjj {I^j{Mjj)) for 
use in Eq. (^. 

The wave function renormalization constants and Z]j^ ^ Eqs. (^ and (^, involve 
similar but more complicated expressions obtained by formally differentiating the expressions 
for n° and with respect to p^, and evaluating the results at = and = M]j, 
respectively. We thus need expressions for all of the self-energy graphs and their derivatives 
with respect to the external momentum, but only for = or = Mfj. 

We have calculated the Feynman integrals using dimensional regularization with the 
number of dimensions D given by -D = 2co' = 4 — 2e. To keep the coupling A dimensionless 
in 2uJ dimensions, we multiply the quartic coupling in Eq. (^ by yU^"^"^ and the trilinear 
coupling by yU^~^'^, where /i is an arbitrary scale mass which will drop out in the final 
results. For the diagrams we need, each loop then appears with a factor which we will 

incorporate into the measure of the momentum integrals, defining Vk as 

Vk = ^4-2- (13) 

Our definitions for the momenta and particle masses in the Feynman integrals are shown in 
Figs. ^ Ordinary Feynman parameter techniques were used to combine propagators, and 



are sufficient for most of the calculations. (However, we used other methods — hyperspherical 
or Gegenbauer polynomial expansions and Kotikov's differential equation method — to eval- 
uate some of the Lemon (C) graphs, as discussed later.) The basic integral needed after 
combining propagators is 



J + 2fc ■ g - m2 + ie]° T{a) [-q^ - m? + it^-^ ^ ^ 

where -M^ + ie = e''^{M'^ - ie). 

Since a number of functions and structures appear repeatedly, we will use some abbre- 
viated notations. Thus, we will scale all masses or momenta by M^, and denote a scaled 
quantity by a caret, e.g., s = s/Mjj. Each loop of a divergent graph contributes a factor 

C = (4V/Mi)^ = {Anfc'Y (15) 

to the final result; we will generally not display this factor except when needed. A derivative 
of an integral with respect to an external momentum, 9(graph)/9g^, will often be written 
as 9(graph). Finally, the first Feynman parameter introduced, x, often appears in the 
combination X = x{l — x). 

Our results for the integrals which appear in the self-energy functions and their deriva- 
tives are given in Appendix A. In the following subsections, we will simply sketch the 
methods used in the calculations, with minimal detail. 

B. The one-loop diagrams 

The calculation of the one-loop graphs is straightforward. However, it is necessary to 
keep terms of 0(e) in the expansion of the resulting expressions for e ^ in order to cover 
cases in which the graphs are multiplied by terms of order 1/e, and for articulated graphs. 

1. Scoop 

The momentum integral corresponding to Fig. |^ can be evaluated directly from Eq. (|14|): 



10 



(16) 



(P - m2) (47r)2 

Q is independent of the incoming momentum, and proportional to the mass in the loop. 
Hence, Goldstone boson Scoops are always zero (m^ = 0) and have simply been ignored 
here. For the case of a Higgs scoop with a self energy insertion, we also need 

„ . dQ ^er(6) 

For the w, Q vanishes identically, and Bq = as well. 

2. Bubbles 

The momentum integral for the bubble diagrams is given in the convention in Fig. ^ by: 

We will need B to order e, an order higher than retained in most calculations. Combining 
the denominators using a Feynman parameter x and performing the T>k integration gives 

^ (Any ' 



a\l-x) + b'x-sX , (19) 
where X = x{l — x) . Similarly, 

a'^{l-x) + b'^x-sX ^ \ (20) 
The cases we will need correspond to the choices 



dq^ {ATr^Mfj Jo 



Bq : s = 1, a = b = 0; 
Bi : d = 1, b = s = 0; 
B2 ■ d = b = s = 1. 

If a = 6 = or if s = 0, the integrals give beta functions and the answers can be written 
down immediately, as can the derivatives dB for the same cases. For B2, we expand the 
denominator in powers of e. The result at overall order e involves the integrals: 
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[ dxln(l - X + x^); I dxhi\l - X + x^)- ( dx ^^^^ x + x^ ^ 

J J J 1 — X + X'^ 

These can be evaluated using the factorization 

l-x + x"^ = {l-c+x){l-c_x), c± = e^^^/l (22) 

The last two integrals then lead to dilogarithms and, in particular, to the imaginary part of 
dilogarithms, Clausen functions |2^. These integrals contribute such characteristic quanti- 
ties as C = a(7r/3) = 1.01494 • • •, vrv^, and TrVSlnS. DeVoto and Duke |2| give a number 



of results which are useful in these and later calculations in an extensive compilation. Further 
results are given in Appendix D. 

We will also need i3o(s) = s^'^Bq and -62(5) for s = ^ Mfj. The second function can 
again be evaluated by expansion in powers of e, or by using a method based on a Barnes'-type 
representation for the denominator in Eq. (pO]). We have found the second method — which 
is apparently new in this context — to be a very powerful way of handling the integrals over 
Feynman parameters encountered in dimensional regularization. The method is discussed 
in detail in Appendix B. 

3. Triangle 

The Triangle graph in Fig. | is related to the Bubble by T = idB/da^, 

^ = /^^((P-a2)) ((fc + g)2_62)- (23) 



Hence, from Eq. (0), 
.dB _ 

'd^~ {An)mjj Jo 



and 



There are four cases of interest here: 
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Tq : s = 1, d = b = 0; 
Ti : b = 1, d = s = 0; 
T2 : d = 1, b = s = 0; 
% : d = b = s = 1. 

The first three lead to beta functions, as for the analogous bubble Bi. By power counting, 
or by comparison to the analogous Bi, it can be seen that the divergences of % and Ti are 
infrared divergences. 7j was evaluated by expanding the integrand in powers of e following 
the method used for B2- We will also need %{s) = s~^~'^Tq and '73(5). The latter was 
evaluated both by the Barnes' method in Appendix B and by direct expansion. 



C. The two-loop diagrams 



1. Sunset 



The Sunset graph in Fig. is one of two new topologies at this order, and is the only 
two-loop self-energy graph in a symmetric 0^ theory other than the double Scoops. The 
momentum integral is given in the conventions of Fig. ^ by 



S J VkVi _ ((k-e + qf - 62) (£2 _ c2) ■ ^^^^ 

Combining the first two denominators using a Feynman parameters x and performing the 
k integration, then combining the remaining denominators using a second parameter y and 
integrating over I gives the expressions 



S = dxX-^ / dyy^-' c'{l - y) + a'y - sy{l ~ y) , (27) 



l-2e 



and 

dS —'if^'^^T(2f] /■! /■! r i-2e 



dq^ (47r)4 Jo Jo 

where = a'^/M]j with 
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a' = - + (29) 
X 1 — X 

The four cases of interest here are 

Sq : a = b = c = s = 0; 

51 : c = s = 1, a = b = 0] 

52 : a = b = 1, c = s = 0; 

53 : a = b = c = s = l. 

The So case is curious because there is no internal scale in the graph: all the masses are 
zero and the integral is to be evaluated at = s = 0. For vanishing internal masses, So is 
proportional to s^"^*^ multiplied by finite integrals, hence vanishes for s — > 0. The case of 
dSo is more delicate since dSo oc s~'^'^ and there is a final infrared divergence unless e < 0. 
We will adopt this limit, with e — > 0^. Alternatively, to avoid this problem, we note that 
So can also be rewritten in a form which contains a mass scale and avoids conflicts between 
the infrared and ultraviolet regularizations by using the identity 

i {k'^ — w?) i 



k? {kP' — m?) k"^ 

i m? 



[k"^ — m?) i (/c^ — m^) k"^ 

This substitution gives 



(30) 



- ^'^''=' -~J ^^^^P (P-m2) (fc^Tp (31) 
with a similar expression for dSo- In either case, the two contributions, calculated indepen- 
dently, are identical, and both So and dSo vanish. 

Since all lines in a Sunset graph are formally equivalent, a judicious choice can be made 
of which to choose as massive in the cases Si and S2- The calculations simplify enormously 
with the choices above. For both Si and ^2, the x and y integrations then factor in Eqs. 
(p7|) and ( pHD and lead to beta functions. 

The all-massive case 1S3 is troublesome, because the last factor in Eq. ( P7| ) mixes singu- 
larities in the x and y integrations. One useful approach, the "(9p" or "partial p" method, 
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makes use of integration by parts in the original momentum-space integral in Eq. (pB]) and 
in the modified integral over Feynman parameters which results to make the divergences 
explicit. It is then possible to expand the remaining integrand in powers of e. The calcula- 



tion proceeds along the general lines sketched for this case by Ramond pOf, but a complete 



calculation including finite parts is what is needed here. The details are given in |2T| . 1S3 can 
also be calculated directly using the Barnes' representation technique discussed in Appendix 
B after a factor {y/ XY~'^^ is extracted from the last factor in Eq. (^). For dS^, the original 
expression in Eq. ( p8|) is sufficient, though the Barnes method is also simple. 

2. Acorn 

As is evident from Fig. the Acorn diagram has the character of a self energy bubble 
with a renormalized coupling. The Aq diagram in Fig. ^ is a Scoop with a self-energy 
insertion. The two types of diagrams give Feynman integrals of the same form with the 
restriction g = for Aq. The common momentum integral is 



A y ^^^^^^^^2_a2) ^^_£)2_52) (£2_c2) ((£ + g)2_^2)- (32) 

Combining the propagators of the Vk loop and integrating, 

^ ^^(47r/i^)T(2-^) .1 ^ _ ( ^ I _^ i , . 

(47r)2 70X2— 7 (£2_^2)2-^ (£2_c2) ((£ + g)2_^2) ^ ^ 

where cP' is defined in Eq. (|29|) . Introducing two more Feynman parameters and performing 
the Dl integration, we obtain the expression 

&y + d\\ -y)- %(1 -y)- (c^ - c^)yz] . (34) 

The particular Acorn graphs which we need correspond to the following sets of parameters 
grouped according to the diagrams in Figs. 2 and 3: 
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A4 : 


a 


= b = 


c = 


d = 


s = 


1; 




A2m'- 


c 


= d = 


s = 


1, 


d = 


b 


= 


Alrn- 


a 


= ^ = 


1, 


b = 


c = 


d 


= 


A, : 


a 


= b = 


c = 


1, 


d = 


s 


= 


A20 : 


h 


= d = 


1, 


d = 


c = 


s 


= 


Aio : 


c 


= 1, 


a = 


b = 


d = 


s 


= 


AiQ-. 


a 


= b = 


c = 


d = 


1, 


s 


= 


A2Q- 


c 


= d = 


1, 


d = 


b = 


s 


= 


Aiq: 


a 


= 1; 


b = 


c = 


d = 


s 


= 



The integral for factors and gives a product of beta functions. The 
integrals can be evaluated by writing the last factor in Eq. as [■ • ■] = [P - 
integrating by parts with respect to z to eliminate the z^^^ singularity: 

1 . . z'Q 



+ 2 / dz- 
Jo 



remaining 

f Qz] and 

(35) 



_e(P + Q)2^ Jo [P + QzY+^\ 

We can now expand the factor z"^ in the last integrand to 0(e) and evaluate the remaining 
integral on z, and find that 



-e^T(2e) dx , 



^ + 2eU2 ( ^] + Oie'^ 



eP2 



(36) 



where Li2(-2) is the dilogarithm function |28], and 



P = c^y + d\l -y)- sy{l - y), Q = -y{& - c?). 



(37) 



In most cases the dilogarithm needs to be integrated by parts with respect to y. This creates 
a set of awkward but finite integrals. 

Differentiating of the expression in Eq. (|36|) with respect to s gives an expression for dA 
with much the same character, 

^2T(1 + 2e) /•! dx /•! 



dA 



(47r)4M 



P \ P 



0(el 



(38) 
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The Aq graphs are Scoops with a particle exchange across the loop; dAq therefore 
vanishes identically. For ^4, P = 1 — y{l — y) and, as for B2, the term in Eq. ( pUf ) 

must be expanded to O(e^), and leads to integrals related to dilogarithms. The same term 
appears in A2m- In all other cases, the term leads to a beta function when integrated, 

and the hardest part of the calculation is the extraction of the finite pieces. 

We will also need Aim, d Aim/ da?, dAim/db"^, dAim/dd?, A2m, dA2m/da? and dA2m./d<P', 
all as functions of s, for our calculation of the Lemon graphs £1, £2™, and their derivatives. 
These integrals were all evaluated using the Barnes' representation method. The results of 
the Acorn calculations are summarized in Appendix A. 



3. Eye-in-the-Sky 



The Eye-in-the-Sky diagrams in Fig. |^ are simply self-energy bubbles with self-energy 

bubble insertions. The most singular terms from these diagrams will cancel exactly against 

the mass counterterm and Scoop insertions. The generic momentum integral is 

/ . \ 2 



S = VkVi- 



and is simply a derivative of an Acorn integral, S = idA/dc^. Thus, from Eq. 
^e'T(l + 2e) 



S 



dxX-' dyy^+' dz z''\l 
Jo Jo Jo 



X c^y + d^{l -y)- sy{l - y) - (c^ - a^)yz 



(39) 



(40) 



Alternatively, 



ie'T{l + 2e) /■! dx /■! ^ 1 , , 



d£ 



{^T^YMl Jo X' Jo 
1 -e 



p 



(41) 



(42) 



17 



when these expressions are well-defined. 

We will need the same six sets of parameters as appeared for the ordinary Acorn diagrams: 





= dc^Ai : 


a 


= b = 


c = 


d 


= s = 


1; 




£3 


= dc2A2m'- 


c 


= d = 


s = 


1, 


d = 


b = 





Si 




a 




1, 


b 




d = 





^4 


= <9c2^3 : 


a 


= b = 


c = 


1, 


d = 


s = 







= 8^2 Alo ■■ 


c 


= 1, 


d = 


b 


= d = 


s = 





S2 


= dc2A2o ■■ 


b 


= d = 


1, 


d 




s = 






These graphs display characteristics similar to the parent Acorn graphs. Thus, S2 = 
idAio/dc^ was calculated using the original Feynman parametrization in Eq. (^Ol) . £^ and 
£^ were calculated by expanding the factor P~2e ^g^j^is in powers of e. No expansion was 
necessary for £2 ^4- ^1 ^^is calculated as a function of s and for s = using the Barnes' 
technique. 

£2 = idA2o was calculated by a different route. The problem here is that the In P/P 
term is ln(l — y)/{l — y), i.e., divergent. The original momentum integral was rewritten 
using partial fractions, as 

((£ + g)2-M^) £2 M]j[{i + qy-M^ M]j {{i + - M^) ' ^ ' 

The last term is proportional to after integration, so vanishes in the limit q^ = needed 
for £2- The first two terms give Acorn integrals, and £2 therefore splits into a difference of 
known terms, 

S2 = J^{A20-A,q). (44) 

d£2 was calculated using the Barnes method. 

4- Lemon 



The Lemon topology shown in Fig. ^ was by far the most troublesome, and three very 
different techniques were used to handle these diagrams |3l[]. In addition, four separate 
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numbers were ultimately extracted by numerical evaluation of two-dimensional integrals 
using the program VEGAS |32]. The momentum integral is 



This integral is ultraviolet convergent, but one case, ^£20 is infrared divergent. The param- 
eter sets needed are: 

£5 : a' = h'' = ^ = d'' = e^ = s = Ml; 
£2^: a2 _ ^2 _ ^ _ ^2^^ ^2 _ ^2 _ g2 _ 0; 

: c^ = s = Ml, a^ = b^ = d' = e^ = 0; 

£20: a^ = e^ = Mfj, = c^ = d^ = s = 0. 
£3 : a2 ^ ^2 ^ ^2 ^ ^2^^ ^2^g2^^^0; 

The momentum integral can be performed formally using Feynman parameters. Com- 
bining the a and b propagators first, adding in the c propagator and integrating over k, then 
combining the result with the remaining d and e propagators in one step and integrating 
over i gives 



C = — ; — -- — 7^ — / dxxX'' dy dudv dw 5(1 — u — v — w) {M — sS) 
{An)^Mfj Jo Jo Jo ^ ^ ^ ^ 



-l-2e 



(46) 



and 



dC = — ; — / dxxX'' I dy I dudv dw u^6(l — u — v — w) 0X9-, (47) 

(AttMh)^ Jo Jo ^ Jo ^ ' {M - sSf^^' ^ ' 

where 



M = u 



dx{l — y) + bxy + c(l — x) + dvX + ewX. 



S = ux'^y{l — y) + X{uy + w){l — uy — w). (48) 

The Feynman parameter integrals are generally too complicated for direct integration 
except for the 11^ cases where s = and just mass terms survive. In particular, £20? <9£2o and 
£3 can all be calculated analytically. (9£3 was calculated partially analytically and partially 
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by a two-dimensional numerical integration. However, the Feynman parameter expressions 
are not generally well-suited to numerical integration. For example, the expression for £3 
had to be reduced analytically to a two dimensional integral before the result obtained 
by numerical integration agreed with the analytic integration. This is presumably because 
VEGAS could not adapt to the complicated singularity structure. 

A more useful technique for £5, the all-massive case, was based on the hyperspheri- 
cal formalism In this technique, a Euclidean propagator is expanded in Gegenbauer 
polynomials (or equivalently, Chebyshev polynomials of the first kind 



1 Zkl 



Y.ZIlCI{K-L), (49) 



{k - if + m2 KL 

where L and L are the magnitudes and four- dimensional unit vectors of the Euclidean 
momenta k, £ and 



m 7^ 0, 



min — , — 1 for m = 0. (50) 



2KL 

fnr 

L' Kj 

The angular integrals are easily carried out for planar graphs using the orthogonality rela- 
tions 



c^^(a-^)c^i(^-£)^ = ;^c^^(a-a). (51) 



Only propagators with angular dependence need be expanded, i.e., for the Lemon, the 
propagators of the 6, c and e lines (see Fig. The result for £5 after performing the 
angular integrations, resumming the series, and making the analytic continuation back to 
Minkowski space, is 

.00 .00 dL^ in(i _ ZkqZklZlq) 



(52) 



-Ml 



g2 Jo (47r)2 Jo (47r)2 (iT^ + a?){L^ + (P) 
There is no known change of variables that will reduce this case to an analytically integrable 
form ||33|. Equation ( ^2]) was integrated numerically instead. The advantage of the hyper- 
spherical formalism is that the integral has been reduced from four to two dimensions, and 
the integrand is a much smoother function. 
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The remaining integrals £i, dCi, £2™, dC2mi and dC^ are all finite, but numerical in- 
tegration of either the Feynman parameter or hyperspherical expressions is complicated by 
singularities or, in the hyperspherical case, a complicated complex analytic structure. For 



these a third technique was used, Kotikov's method of differential equations Kotikov 
uses the method to find a first order differential equation in the particle mass with, hope- 
fully, "easily calculable" inhomogeneous terms. In the present case, we obtain a differential 
equation in s (the only free variable is the ratio s = s/Mjj) by taking the derivative of the 
momentum integral in Eq. (46) with respect to the incoming momentum. The equation can 
be simplified by using the six identities 



dkj di^[p^{---)] = Q (53) 

where (■ ■ ■) is a product of propagators, and is any of the variables /c^, and g^, and 
is either or C^. 

The result is an inhomogeneous first order differential equation with the inhomogeneous 
term expressible in terms of Bubbles^ Triangles^ Acorns, and their derivatives. C can then 
be found by direct integration. This is an interesting technique, and the remaining cases 
were handled this way as discussed in Appendix C. Our exposition there is less general than 



Kotikov's 34 since we treat only specific diagrams. 



D. The self-energy terms and renormalization constants 

We will write the expansions for the self-energy terms and wave function renormalization 
constants in powers of the coupling A as 

= n(i) + n(2) + . . . 

Z = l + + + . . . (54) 

The one-loop self energies are given by the first two terms in Eqs. (|TT]) and (|T^ with the 
wave function and coupling renormalizations ignored. 
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- + (1 - ^) + (2 - 27 + 7' + C(2)) + 0(e^ 



(55) 



and 



nff(Mi) 



3A 



+ (9-71^3-57) 



+e (^^7^ - 97 + 17 - y - 2nV3{l - ^7) + TrVsinS - AVscj 

+27r + (2 - 7)i7re + 0(e2) . (56) 

Here 7 = 0.5772- ■ ■ is Eulers constant, C(2) = n^/G, and C = Ci (f) = 1.01494- ■ -. We 
have not expanded the factors C,'^ for simphcity. This will also make the disappearance of 
the arbitrary scale fi manifest in our final results: the complete cancellation of the poles in 
e will allow us to take the limit e ^ 0, ^'^ = {4tt fi"^ / MfjY 1 with no pieces left over. The 
one-loop contributions to the wave function renormalization constants are similarly given by 

A 



167r2 



r l + e(|-7) + 0(e^ 



(57) 



and 



.(1) 



167r2 



r fl2 - 2ttV3 + e(24 - I27 - 371^3(1 - I7) + 27ry3 In 3 - 8^30) + 0(6^ 



(58) 



The leading terms in these expressions agree with the results of previous one-loop calculations 

The real parts of the self energies at O(A^) involve the corrections to the one- loop self 
energies from the one-loop wave function and coupling renormalizations as well as the new 
two-loop contributions. The algebraic programming system REDUCE was used to manipu- 
late the expressions in Eqs. (0) and (|12D and obtain the expansions needed. The results, 
evaluated on mass shell, are 
A2 



(IGtt 



2^2 



e^M^H ( + -(-22 + 157 + 3v^7r) 



-72^5 In 



V5 + l\ ..,^2 /^\/5 + l\ , ool^o on 2 



24 In^ 



32 In 2 - 9O7' 
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and 



753 

+2647 + 40C(2) - ^ + 66 VSC - 36V37r7 + 45^371 
-187rv/31n3 + 0(e) ). (59) 



Reng)(M^) = -^e'^M^f^ + -(-71 + 427 + lOv^: 



-1441n2 - 457r2 -3787^ 

-367^2 - 162ir5 + 12787 - 63C(2) + 90C(3) 

-1524 + 126\/3C - 180y37r7 + 408V37r 

-81^(2)ln2 - 907rV31n3 + O(e)y (60) 



The constants Ki are from numerical integrations, and are given in Appendix A. These 
numbers were calculated to accuracies which vary from a part in 10^ to a part in 10^, more 
than sufficient for our purposes despite the large topological weights and the cancellations 
which occur. 

The O(A^) contributions to the wave function renormalization constants obtained from 
the derivatives of Eqs. (|TT|) and (|I^) using REDUCE are 



(167r2)2^ \ t V 2 ; V 2 

100 38 , , 152 

— n 2 + 67 C 2 - 12/^3 

9 ^ 3 ^ 54 



In 2 + 67 - ?C(2) - 12/^3 - ^ + 3v^7r + 0(6)^ (61) 



and 



^^'^ = (1^^'' ("7 + 144 In 2 + 427r2 + 81C(2) In 2 

+ 36^2 + 162^5 + 67 + 33C(2) - 90C(3) 

273 \ 
-16^371 In 3 + 334 VSC + — - 292^371 + 0(e) j . (62) 

The leading 1/e divergences in Z^^^ and z'^^ are identical because of the original S0(4) 
symmetry of the theory which is preserved in the interactions at large momenta. 

Finally, the bare coupling Aq can be identified as the coefficient of — |(wq + H^)'^ in Eq. 

(0), 
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Ao = A(i_5£5M^gWaMy (63) 
Substitution of the results above using REDUCE gives 

Ao = A + (— + 25- 127 - sVSn + e [37^^/31n3 - 12v^C - 67rv^(l - ^7) 

- vr^ + 6f - 257 + f ] + 0(e^)) + ^ + ^(29 - I67 - 4^3.) 
32906 

+ I2K3 + 162/^5 + 36/^2 - 90C(3) 
+^C(2) - 10447 + 2887^ + 547r^ + In 2 

+81C(2)ln2 - 3697r + 1447r\/37 - 60V3C 

+727r V3 In 3 + 0(e)^ + O(A^) . (64) 



IV. THE SCATTERING AMPLITUDES 

The two-particle scattering diagrams generated by the quartic interactions in Eq. (7) are 
shown in Fig. 4. There is just one generic one-loop diagram with an associated amplitude 
Bs, and two two- loop diagrams with amplitudes Vs and ^5, the last named after the similar 
Acorn self-energy diagram. The small number of diagrams is the result of the suppression 
of the dimension-three trilinear couplings in the limit ^/s 3> M^, as was noted earlier. The 
functions Bs, Vs, and ^5 are given in Appendix A. 

We note that Mh can be neglected on internal hues in the high energy hmit, and that the 
scattering amplitudes consequently have the SO (4) symmetry of the bare quartic couplings. 
This symmetry is reflected in the structure of the 4x4 matrix J- of unrenormalized transition 
amplitudes. All of these amplitudes can be expressed in terms of a single function A{s, t, u), 
where 

A{s, t, u) = -2Ao - i(-iAo)'A« - i{-i\ofA^^\ (65) 
A^'\s,t,u) = 16Bs{s)+iBs{t)+iBs{u), (66) 
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A^^\s, t, u) = lOAVs(s) + 8Vs{t) + 8Vs(u) 

+176Asis) + 80^5(^) + ^OAs{u), 



(67) 



and Ao is the bare coupling defined in Eqs. (63) and (64). We will write the two-body 
channels in the order w~^w~, zz, HH, and zH. then has the structure 



V 



(68) 



^(s) + ^(t) A{s) A{_s) 

A{s) A{s) + A{t) + A{u) A{s) 

A{s) A{s) A{s) ^ A{t) ^ A{u) 

A{t) 

where we have indicated only the first variable in A since this function is unchanged by an 
interchange of the remaining two variables. 

As noted after Eq. (62), the leading terms in 1/e in the renormalization constants Z-^ 
and Zh are identical because of the SO (4) symmetry of the initial theory. We will therefore 
make an intermediate renormalization which preserves the SO (4) structure by multiplying 
A by the field renormalization factor {Z]j^Y replacing Aq in Eq. (65) by the expression 
in Eq. (63). This gives 



AR{s,t,u)^-2\\\- 



Ren2,(M^)-no(o)- 



Ml 



+^A^^« {^_^^^Sm^&^^y^^A^^,^ (69) 

where each term is to be expanded to 0(A'^). Ar is just the physical amplitude needed to 
describe the scattering in the w+w", zz sector. Channels involving external Higgs bosons 
are multiplied by a finite renormalization factor {ZhIZ^^'^ for each external Higgs boson. 

After a number of nontrivial cancellations, A^ reduces in the limit e ^ to a finite 
expression independent of the arbitrary scale which appeared at intermediate steps of the 
dimensional regularization, 

\2 

AjAs, t, u) = -2A + -— f-16 In(-s) - 4 Ini-i) - 4 In(-'a) + 2 + 6v^7r) 
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A3 



(167r2)2 



^ -1921n2(-s) + 1761n(-s) + 96v^7rln(-s) ^ 
-481n^(-t) + 801n(-t) + 24^371 In(-t) 
-481n^(-a) + 801n(-ji) + 2A^/^n\n{-u) 
+60a/5 In ^ - 456 V3C + 138A/37r 
+240 In^ ^ - In 2 - 1807r2 



V 



(70) 



-72^2 - 1^C(2) + 180C(3) 
-324^5 + 24^3 + ^ - 162C(2) In 2 

where —s — e^^'^s, and the quantities —t and —u are real and positive in the physical region. 
Substitution of Ar for A in Eq. (69) gives the partially renormalized transition matrix 
the fully renormalized matrix of (Feynman) transition amplitudes is 



M = ZJFoZ 



(71) 



where Z is a finite diagonal matrix of ratios of renormalization constants, 



Z = dmg{l,l, Zh/z^,z]{'/zI/'), 



(72) 



and the products in A4 are to be expanded to 0(A3). The SO (4) symmetry of is lost 
in Ad. 



V. COMMENTS 

The principal results of this paper are the expressions wc have obtained for the renor- 
mahzation constants Zyj and Zh, the self-energy functions n° (0) and Ren^(M^), the bare 
coupling Aq, and the renormalized two-body scattering amplitude Aii{s,t,u) at two loops. 
We will use these results in a following paper to study limits on the range of validity of 
low-order perturbation theory in the Higgs sector of the standard model. 

More generally, the quantities connected with the on-mass-shell renormalization of the 
theory can be used in other calculations of W^,Zl scattering based on the use of the 
equivalence theorem, for example, in possible calculations to extend the analysis of the 
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scattering in the "low-energy" limit <^ ^ M}j to two loops. This extension 
would be of considerable interest with respect to studies of the energy at which the effects 
of a very massive Higgs boson would become evident in future experiments p6| . 

We note finally that the new Barnes'-type method we have introduced for factoring and 
evaluating integrals over Feynman parameters has a wide range of applicability as illustrated 
in the calculations above and in Appendix B. The results collected in Appendix A and 
Appendix D should also be useful in future calculations. 
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APPENDIX A: RESULTS FOR THE SELF-ENERGY GRAPHS 

This appendix is a listing of the results of the Feynman graph calculations. Some in- 
termediate integrals used in the differential equation method are also listed. Since only the 
real parts of the self-energy graphs were needed for the counterterms, strict attention was 
only paid to the real parts of those integrals when expansions were made. 

A factor of ^YlSvr^ should be included for each loop in a graph. Factors of Mfj should 
be included so that, with the couplings given in Eqs. ( pH]) and (p!2D, the self-energy graphs 
have dimensions of mass squared. That is, Q, 5 oc M^; B, ^ oc 1; and T, £, C oc Mjj^ . 
The derivatives of these functions will have one power of M'jj less. Note that the graphs 
with a Q subscript have no dependence on the external momentum, so c}(graphQ) = 0. 

Expressions which involve explicit expansions in powers of e have been evaluated using 
expansions in e in intermediate steps. These expressions are not correct beyond the highest 
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power of e indicated. Terms involving gamma or beta functions have been kept in the 
forms given for compactness, and are to be expanded in powers of e when used. These 
terms can frequently be rearranged, and may be replaced in some cases by quite different 
expressions which are equivalent to the relevant order in e if an intermediate expansion is 
done differently; only the final expanded expression is meaningful. 
In the following, B{a, h) is the beta function, 

T{a)T{h) 



B{a,b) 



T{a + h) 



Li„(x) is the polylogarithm [2S 

rx 

Un{x) = / Li„_i(t) — 

JO t 

and ({n) is the Riemann zeta function. 
Constants: 



oo 



71' 



(Al) 



(A2) 
(A3) 



Scattering graphs: 



7 = 0.57721566... 

C = Ci{n/3) = 1.014942 
K2 = -0.86518 
K'^ = 0.548311 
K3 = -0.1066639 
K5 = 0.9236306 



Bs = -t{-s)-r{e)B{l-e,l~e), 



Vs = {Bs)\ 

As = -{-s)-^T{2e)B{2 - 4e, e)B{l - 2e, 1 - 2e)B{l - e, 1 - e) 



(A4) 
(A5) 
(A6) 



Scoops: 
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Bq = -ir(6). 



Bubbles: 

Bo{s) 
Bo 

B,{s) 
dB^(s) 
B, 
dBi 

Us) 

B2 
dB2 
Triangles: 

Us) 
% 

dr, 

% 

dT2 



-z(-s)-r(e)S(l-e,l-e), 
-ie^^T(e)5(l - e, 1 - e) 

-^r(e) [1 + e (^2 + i(l - s) ln(l - s)) + O(e') 
-ir(l + e) ln(l - - ^ + 0(e)) 

. r(i + e) 



'(l-e)(2-6)' 



4—5 s 

-ir(e) ( 1 — 2e , = arctan 



\/4s - s 



' V4s - s 

+ ]-e^ [ dx In^ [1 - sx{l - x)] + 0{e^)) , 
2 Jo J 

-iV{e) 



+ 2e 



/ TT \ ^1 27r TT In 3 4C 



ir(l + e) 



27r 



e 2-^ 



V3 

TT 27rln3 8C 



V3 3^3 3^3, 



+ 0(6^) . 



_(_5)-i-T(l + e)5(l-e, -e), 
e^''T(l + e)B(l - e, -e) 
-e^'^T(2 + e)S(-e,l-e), 

r(6) 

1-e' 

-r(2 + e)S(3,-e), 

_r(l + e) 
1-e ' 

r(2 + 6) 

(l-6)(2-e)' 
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r3(s) = -r(i + e) 



arctan 



^ xln[l — sx{l — x)] 



£(1 + 6) 

3^3 
r(2 + e) 



[7r-e(7rln3-4C) + 0(e2)], 

Un 7rln3 4C 



TT 



3^3 



+ e 1- 



+ 



3^3 3^3 3^3 



+ 0(6^ 



Sunset: 



(A26) 
(A27) 
(A28) 



50 = dSo = 0, 

51 = 2i^^5(l - e, 1 - e)5(l + e, 2 - 4e), 



95i = -2*^1' 



^2 = I 



dSo = -i 



(l-26)(l-6)' 

.r(2e)E(l + e,l + e) 



(2-e)(l-6) ' 



..[r(6)f 



l-e VI -2e 
1 r^,. .,2 .3 



(l + 26)(2-6) 



953 = -i-[r(l + e)]^ + i- + 0(e), 



>1 corns; 



= -r(2e)S(e, 1 - 2e) 

e 



5(1 -e, 1 -e) 



5^20 



2-e 



10, 



r(6)r(26)r(i-6) [r(.)r(i-6)f 

1 -e 
r(e)r(l + e) 



-^ + 0(e), 



l-e 
r(2e) 



(1 - e)r(2 - 2e) 
S(l-e,l-e) + 2-C(2) + 0(e), 



e''^^ + 2C(2)-3 + 0(e), 



5(1 - e, 1 - e)5(l + e, 1 - 2e) + 2C(2) + 0(e), 



r(l + 2e)5(l -e,l -e) 



(2 + e)e 
Am = A-2C(2) + 0(e), 



-2 + C(2) + 0(e), 



1 2C 



(A29) 
(A30) 
(A31) 

(A32) 
(A33) 
(A34) 
(A35) 



(A36) 
(A37) 

(A38) 
(A39) 
(A40) 
(A41) 
(A42) 
(A43) 
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^3 = 

A = -r(2e) 



dAio 



V{2e)B{l - e, 1 - e) 



(l-e)6 
r(l + 2e)5(l -e,l -e) 



(2-e)(l-e)e 
5(1 -e,l -e) 



+ 2A/3C + 0(e), 
+ l + 0(e), 



TT 



, . +2S(l-e,l-e) 2 ^ 



/ 27r 7rln3 1 7C\ , 



9A = 



r(l + 2e)5(l-e,l-e) 



27r ^ 



/ TT 27rln3 1 8C \ , 



r(26) 

e(l-6 
r(2e)S(l-e, 1-e 



r(i + e)r(i-6), 

+ l-C(2) + 0(e), 



6(1 -e) 

r(2e)5(l-e, 1-e) 



0, 



^ = 1,2,4. 



2C 

+ l + -^ + 0(e). 



Auxiliary Acorns: 

Am(^) = r(e)[r(l-e)] 



2/ r(2e) 



r(6)(-^)- 



dA 



Im 



dd? 



r(2 - e) (1 - e)r(2 - 2e) ; 
+ (l + J) Li2(s) + 2 (^1 - ^) ln(l - s) - 3 + 0(e), 

- In(-s) + (l - t) ln(l - s) + jLi2(s) + 0(e), 

[r(e)r(i-e)f 



(1 - e)r(2 - 


-2e)^ 


+r(-e)r(i - 


oo 
n=0 


[r(e)r(i- 




(i-e)r(i- 


-2e)^ 


-r(-e)r(i - 


oo 


[r(e)r(i- 




(l-e)r(l- 




-r(-e)r(i - 





(n + 1 + e)r(n + e)r(n + 1 + 2e) 

r(n + 2)r(n + 2-e) * ' 

r(n + l + e)r(n+l + 2e) 
(n + 2)r(n + 2 -e)r(n + 1)^ 



-s) 



+ TLi2(s) 

5 
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1 

+- 

s 



1 + ( 1 - ^) ln(l - s) - 2 (^1 - i ) ln'(l - s) 



+ 0(e), 



^•^2m(^) = r(l + 2e)S(-e,l-e) 



1 - ^ ) ln(l - s) 



2(4 - s) 
H — , _ arctan 



VAh - s 



S2 



V4s - 



- y" (^ln^(A - ^/) - InM + Iny [In A - ln{A -y)]- Li2 
+0(e), where A^l- %(1 - y), s < 4, 



^•^^-(5) = lr(l + 2e)i?(l-e,l-6)- ^ 



arctan 



[ dy 

Jo 



+ 1 dy — r^[lny -\nA-\n{A-y)] + 0{e) 
BAa . . 27r 7rln3 7C 



■in-the-Sky: 



= r(2e)S(l - e, 1 - e) 



[r(e)r(i -e)]2 , 1 , . /I, / 



ln(l - s) 



dSi^-i{l + e) 



(l-e)r(l-2e) 

[r(6)r(i-e)p 



2s s 
1 



::Li2(s) +0(e), 
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(i-6)r(i-26) 

/r(2 + e)r(e) 



2-C(2)-^-j+0(e), 
e- + i[C(2)-l] + 0(e), 



5(3,-6) + 2C(2)--J+0(e), 
ir(l + 2e) ^^^7''^~'^ ^(e, 1 - 2e), 



1-e 



9^2 



1 + 2e 
2(2^ 



^^3 = <?5 + ^7C(2) + 0(e), 




1 I 



d£z = ir(2 + 2e)B(l - e, 1 - e)- - 



TT 



TT^ 7rln3 7C tt ' 



+ — - 



+ 



3 12 9^3 9^3 V3 



+ 0(6), 



^^^T(l + 2ey(l-e,l-e)^ /^_4C\ ^^^^^ 



(l-e)e 



V3, 
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Lemons: 



r(2 + 26)i?(l-6,l-e) .3 



2% 

TT In 3 7C TT 



£5 = ir{2e)B{l - e, 1 - e 



--C(2) +0(e), 



3v^ 3^3 VS 6 
a£5 = ir(2 + 2e)B(l-e,l-e)- ^ 



TT 



e V3 9a/3 
2 7rln3 31C tt ^ 



3 36 9y/3 9^3 v^, 



+ 0(6). 



dC3 



2Li3(-s) - 2Li3(s) - 4Li3(l - s) - 4Li3(l + s) 
2 liisLi2(-s) - In^ slii(l + s) - 7rMn(l + s) + C(3) + 2i7rln^(l + s) 



+41n(l-s) (Li2{l + s) + -C{2) 

>Ci = -i6C(3) + 0(e), 

9A = ^[6C(3)-2C(2)] + 0(6), 

£2o = i[C(2) + 41n2] + 0(e), 

r(2 + 2e)B(l + e, 1 - e)S(3, -2e) 



+ 0(6), 



(2-6)(l-e) 

z (4 In 2 - ^C(3) - Im + Im In2 + K, + 2mK'^ + 0(6) 



/ , , 2C 27r2 
-/:2m + MC(2)-^ + i 



-41n2-C(2) + 



V3 3^3^ 
6V51n 



0(6), 



+ 2 In' 



0(6) 



4-41n2- -C(2) + 81n2 

2 \ ^ 

~- In — + + 0(6) 
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dC^ = -i (^7^5 + ^C(2) + ^)+ 0(e). 



Auxiliary Lemon integrals: 

K2 = 



K' 



1 41ns 

as , arctan 



V4 



s — s^ 



+ ds dy- 

JO Jo 2 — s 



s — s" 



^ Iny ln(l — sy) 



— ln^(a — y) — i In^ a + In y In 
s V 

a = 1 - sy{l - y), 



a 
a 



Li2 - 



S 2 s 
as = arctan ■ 



a/4s - 

"dx^dy^'^'-y^ 



V^s^ 



Jo 1 — xy 



1 



^ 2 x^{l-xYil-yy 
4:x{l - x) + x^y + - ^ ) y 



xln 1 + 4^ -3- 
Y x\ \ — x) ) 



1 - xy 

x{\ - x){\ - y) 



1 — xy 



1 — xy 



/■I /■! 1 / b b 
K,= dx dy — In {1 + ^(1 + a, + 
Jo Jo xy \ 2 



(A82) 



(A83) 



(A84) 



(ASS) 



% ^ V (1 + (^x + dyY - 4:a^ay) , 



(A86) 



APPENDIX B: USE OF A BARNES'-TYPE REPRESENTATION TO EVALUATE 
INTEGRALS OVER FEYNMAN PARAMETERS 

We have found it very useful in evaluating a number of the integrals over Feynman 



parameters to use a Barnes' type representation for the binomial expansion [ p7| 

1 rioo 

T{a)il- = — dtT{t + a)T{-t){-zy, |arg(-2)| < vr, (Bl) 

2711 J ~ioo 

where the contour of integration runs to the right of the poles of r(t + a), i.e., t = —a, —a — 
1, . . ., and to the left of the poles of r(— t), i.e., t = 0, 1, . . .. The binomial expansion of 
(1 — z)^"" for |2;| < 1 (or \z\ > 1) is recovered when the contour of integration is closed at 
infinity in the right half (left half) of the complex t plane. The advantage of the integral 
representation is that it holds for any magnitude of and any a. 
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We will illustrate the use of the Barnes' representation by example. We first consider 
which is given up to a factor (47r)*^ by the expression in Eq. (34) with a — 1, b — 
c = d = 0: 

-2e 



yz 

-sy{l-y) + — 

X 



Am(s) = -r(2e) / dxx~\l-x)-' f dyy' f dz z'' 
Jo Jo Jo 

^ -r(2e) j^dxx^l- x)-' £ dy y-' dz z'^'^ [\ - ^"^^^ ~ j . (B2) 
We can use the Barnes' representation for the last factor to re express A\m{s) as 

\ rioo 

Aim{s) = -— / dtr{t + 2e)r{-t){-sy 

X r dxx^+'{l-x)-' ^ dyy-'{l-yY ^ dz z'^-'-^ (B3) 

^0 ^0 ^0 

The integrals on x and y converge and can be expressed in terms of F functions for Re i > 
— 1, e — > 0+. The final integral converges for Ret < —e. We find that 

A (^.) [r(i-6)P p r(t + 26)r(t + 6)r(-t) 

where the contour of integration must cross the real axis between —1 and — e. If we take 
|s| < 1, we can close the contour of integration to the right and find that 

A -T(e)\T(l e)]'(-^^ r(6)(-g)- \ 
A,m-T{e)[T{l-e)] l^f(2^ " (1 - 6)r(2 - 2e) ) 

+ ^,rin + 2-enn + 2f " ^^'^ 
The series is nonsingular for e — > and can be summed in that limit to obtain the expression 
in Eq. (A52), or treated more generally in terms of hypergeometric functions for e 7^ 0. 

This calculation illustrates a fairly typical result, that we can use the Barnes' repre- 
sentation to at least partially factor a multidimensional integral over Feynman parameters, 
and reduce the calculation to one of contour integration. It was essential that the repre- 
sentation held in this case for arbitrary values of sx{l — y)/z subject only to the constraint 
|arg(— s)| < TT. This requires that —s — e~*'^s for s identified as usual with s -\- ie. 

In some cases, partial factorization of the integrals is sufficient. For example, if we write 
the expression for dAimlda^ as 
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•^^""(s) = r(l + 2e) dxx-Hl-x)~' dyyHl-y) dz z'^^' 
Jo Jo Jo 



-l-2e 



X[l-%(1-.)1— (l-^-^j (B6) 



and apply the Barnes' representation to the last factor and integrate over z, we obtain a 
simple series in the variable y/ [1 — sy{l — y)]. The singular term in e can be extracted and 
the remainder summed in the limit e ^ 0. The result is the expression in Eq. (A57). 

The method can also be used to advantage to make multiple expansions. For example, 
if the propagators in the expression for Si{s) arc combined in an order different from that 
which leads to the expressions in Eq. (34), one obtains the awkward expression 



=ir(l + 2e) dxx^+Hl-x)-' j^dyil-y) Cdzz^il-zY 
Jo Jo Jo 



I v-L - V) I y-L- 
Jo Jo 

-l-2e 



X (l - £^^i(lzMj . (B7) 

After applying the Barnes' representation to the last factor, one obtains an expression which 
involves (1— yz)*. This may be factored by a second application of the Barnes' representation. 
After the x,y, and z integrations are performed, the contour integrals are easily evaluated, 
and lead to the result in Eq. (A59). A simpler calculation based on Eq. (40) gives the same 
result. 

Finally, it is not always necessary or even desirable to evaluate all the integrals at an 
early stage even if that is possible. An example is provided by 882 which is given by the 
integral 

d£2^ir{2 + 2e) dxx-'{l - x)-' C dyy'^+'{l - y)'^-^' dz z'^+'il - z) 
Jo Jo Jo 

/ N -2-2e 

Use of the Barnes' representation on the last term gives a complete factorization of the 
integrals over the Feynman parameters. However, the singularities at y = 1 in the resulting 
expression cause difficulties if the y integration is performed at this stage since convergence 
of the y and z integrations require different signs for Ret or e. However, it is clear from 
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the original expression that the apparent singularities in y are absent. In this case, we 
therefore integrate only over x and z. After closing the contour of the t integration in the 
left half plane, we obtain a series in positive powers of (1 — y). The singular term in e can be 
extracted, the remainder of the series summed to simple functions for e = 0, and the final 
integration on y performed to obtain the result in Eq. (A63). 

APPENDIX C: AN APPLICATION OF KOTIKOV'S METHOD 
TO THE LEMON GRAPHS 

In this Appendix we describe the application of Kotikov's differential equation method 
3JI to the Lemon graphs which appear in . The method is based on the identity 

= t'l VkVi—{p^AaA,A,A,A,), (CI) 

where the propagators in the Lemon graph will be denoted by iAa, iA^, ... in the conventions 
of Fig. 1^. The notation C[ ■ ] will be used to indicate the action of the Lemon operator in 
Eq. (46) on the quantity in the bracket, i.e., multiplication by the Lemon integrand and 
integration over the loop momenta k and i. For example Ci[l] = Ci while £i[l/A(j] = iAim- 
If we choose = in Eq. (CI) and perform the derivatives we obtain the identity 

= £ [2cj - 2P A„ -2k-{k + q)Ai,-2k-{k- £) A J , (C2) 

where 2u = A — 2e. The scalar products can be expressed in terms of inverse propagators 
and masses, e.g., 

2k-q = ^-^-q'-a' + b'. (C3) 



The identity can then be rewritten as 



(4 - 2uj)C = C 

-A,{~^^+q'-a'-b^)-2a'A, 
If we choose p^ = in Eq. (CI), we obtain instead 



(C4) 
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= £ 



A.(i,-i, + c^ + rf^ 



-A, 



1 



1 

A, 



+ — 



(C5) 



Further identities can be derived by replacing d/dk^ by d/di^ in Eq. (CI). These amount 
to interchanges of labels in the relations above, with (P and 6^ <-> e^. Several of the 

Lemon graphs are symmetric under this interchange, or under b"^ corresponding to the 

interchange of momenta k — > ~k — q,i —i — q which leads to further identities. 

The result of differentiating a Lemon with respect to the square of the incoming momen- 
tum is the identity 



a dC 1 ^ 

- —^C \-2q ■ {k + g) Ab - 2g ■ (£ + g) A^ 



dq^ 2g2 dq^j, 2q 



This can be rewritten as the differential equation 



dC C 



dq' 



+ 



2g2 



A, 



(C6) 



(C7) 



A, ^ " ' Ae 

The key observation at this point is that the terms in the expressions above which involve 
ratios of A's are related to derivatives of Acorns or Double bubbles. For example, the ratio 
Aa/ Ah on the right hand side of Eq. (C7) removes the b line in a Lemon, thus contracting two 
vertices, and introduces an extra a propagator. The result is C{Aa/ A^) = idA{a, c, d, e)/da? , 
i.e., a ^ a,b ^ c,c ^ d, and c? — e in the Acorn graph in Fig. |^. Terms of the form 
(mass)^ X A are derivatives of Lemons. The problem is to use the identities above and 
the symmetries to eliminate the unwanted Lemon terms in Eq. (C7) to obtain a differential 
equation for C with the inhomogeneous term expressed in terms of known graphs. The 
results for the cases of interest are: 



ds s 

dC2m _|_ ^2m 

ds S 



■A, 


A,,- 




A, 


A,. 


4( 



2-s 
1 



^2m 



. dAA 

'dd^J' 

' ' Ad A. A, Ae. 



(C8) 



1 / A, Ad 
Afc ^ A„ 

dA2m 



dAi 



+ 



m2(2 



-{T,Bo + B2%, 



^ dA2m ^ dAir. 



dd^ 



(C9) 
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and 



1 ^ r A, A AJ 

— ^■'-'1 1 r 




ds s 



1 + s HsAft Ad A, 
i ( 1 dAim _ dAim \ 



l + s\s dh^ dd? I 




(CIO) 



All the terms on the left are finite for e — > 0, but most of those on the right have single or 
double poles. 

The differential equations have general solutions of the form 



where f{s) is the inhomogeneous term in the differential equation. A careful study of the 
limiting behavior of the integrands for s ^ with fixed e > shows that sC{s) — > for 
s ^ for the graphs which appear in Hh . (We did not use this method for the graphs 
which appear in II^^^ which must be evaluated for s — 0.) We can therefore take the lower 
hmit of integration in Eq. (Cll) as zero, and write the £'s as 



where an overall factor of l/M^ has been suppressed following the convention of Appendix 



The homogeneous terms need to be calculated for arbitrary s. This is trivial for Bo{s) 
and Tq{s). B2(s),7^(s), and the derivatives of the Acorn graphs with respect to internal 
masses were all evaluated using the method based on Barnes' representation for the binomial 
expansion discussed in Appendix B and illustrated there for Aim{s), and also by direct 
expansions in powers of e using the expression in Eqs. (19), (24), and (36). The singularities 
in the individual terms in the /'s cancel explicitly, and the problem was reduced to one of 
integration. 

Some integrals in and djC2m were evaluated analytically. The rest, while perhaps not 
insuperable, were evaluated numerically. The quantities actually needed are = ^2m{^) 
and dC2m — {dC2m/ds){l). The second can be evaluated immediately using Eq. (C12), 




(Cll) 




(C12) 



A. 
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dL 

dC^—{l) = -C{l) + f{l). 



(C13) 



The results are given in Eqs. (A77) and (A78) 



Ci was considerably easier to calculate, and was evaluated analytically. Curiously, 



£1 = £0 = -^6C(3) 



(C14) 



where Cn is the result for all internal masses zero and s = 1 |09i, but 



dCi = dCo-t2C{2). 



(C15) 



Finally, dC^ was evaluated using Eqs. (C8) and (C13) and the result for £5 obtained using 
the hyperspherical method, Eq. (52). 

APPENDIX D: A COLLECTION OF USEFUL INTEGRALS 



DeVoto and Duke ||29[ give a useful compilation of integrals involving logarithms and 



polylogarithms which occur in integrals over Feynman parameters. Lewin |2^ is the standard 
reference on polylogarithms. We give here a further collection of integrals which occured 
repeatedly in the calculations of the finite terms in the self-energy functions. The first few 
are elementary; the remainder are not. 
Integrals involving X = x{l — x): 

1 271 



1 1 _ 2 Air 

'^''(1-X)2 -3 + 973' 



dx In(l-X) = -2 + ^, 
^0 a/3 

f\ , 2n ^^ o 47r 27rln3 8C 

/ dx ln^(l - X) = 8 ^ 

Jo V3 

/ dxx ln(l — X) = — 1= ■ 
Jo ^ ' 2^3 



V3 V3 VS' 



(Dl) 
(D2) 
(D3) 
(D4) 
(D5) 
(D6) 
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dx Inx ln(l - X) = 4 - ^ - VSC - ^ 



/ dx 

Jo 



2 

V3 36' 
i^^ln(l -X) _ 27rln3 8C 
^ 1-X ~ V3 ~373' 
1 , In(l-X) -7r2 
X 18 
1 / In(l-X) ^ 1 \ _ _^ ^ JJL 
\ X/ a/3' 

1 Inx _ -2C 

1 xlnx — C TT^ 

ax = — ^ H , 

1-X 36' 

1 Inx -TT 4C 

dx ■ 















(1-X)2 3^3 3v^' 
1 xlnx — 2C 
(1 - A')' ° 371 ' 

Note that many further useful identities can be obtained by substituting 
also by using the identity 



C dxxf{X) = \ f'dxfiX). 
Jo 2 Jo 



Integrals involving the dilogarithm ||2^ Li2(-2): 



t dxU2{x) = C(2) - 1, 
Jo 

I dx lnxLi2(x) = 3 - 2C(2), 
Jo 

I cixln(l-x)Li2(x) = 3-2C(3)-C(2), 
Jo 

j'^dxU^i^^ =-ic(2)-21n2, 



1 27r 
dxLi2(X) = h ^ - 4, 



' 18 V3 

/;cixLi2(i£^)=-2v/3C 

1 /Li2(X) 1\ TT 



dx 



V ^/ 73 9 
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FIGURES 

FIG. 1. The one-loop self-energy graphs needed in the calculation of 11^^ and 11^^ The heavy 
lines represent Higgs bosons, and the thin lines, Goldstone bosons. Tadpole graphs cancel exactly, 
and are omitted. 

FIG. 2. The two-loop self-energy graphs needed in the calculation of 11^^ The weights with 
which individual graphs appear in 11° can be read off from Eq. (11). 

(2) 

FIG. 3. The two-loop self-energy graphs needed in the calculation of II ■ The weights with 
which individual graphs appear in 11^ can be read off from Eq. (12). 

FIG. 4. The scattering graphs through two loops which involve only the dimension-four quartic 
couplings and give finite or logarithmically divergent contributions to the renormalized scattering 
amplitude for s, —t, —u — > oo. Particle masses are irrelevant in this limit so all internal lines may 
be taken as massless. 

FIG. 5. Examples of scattering graphs which contain dimension-three cubic couplings and give 
vanishing contributions to the renormalized scattering amplitude for s, —t, —u — > oo. For example, 
the scattering Eye graph £i{s) with the topology of (d) given in Eq. (A59) vanishes as \v? s for 
s — > oo after renormalization. 

FIG. 6. The conventions used for the particle masses and momenta in the one-loop self-energy 
graphs. 

FIG. 7. The conventions used for the particle masses and momenta in the two-loop self-energy 
graphs. 
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